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Abstract. We discuss the possibility of preparing highly entangled states by simply 
cooling atoms into the ground state of an applied interaction Hamiltonian. As in 
laser sideband cooling, we take advantage of a relatively large detuning of the desired 
state, while all other qubit states experience resonant laser driving. Once spontaneous 
emission from excited atomic states prepares the system in its ground state, it remains 
there with a very high fidelity for a wide range of experimental parameters and 
£h ' all possible initial states. After presenting the general theory, we discuss concrete 

applications with one and two qubits. 

■ 1. Introduction 

Dissipative systems are dynamical systems which lose energy over time. In classical 
physics, this typically happens due to friction or turbulence. In quantum mechanics, 
the loss of energy is due to processes like the spontaneous emission of photons. The 
Q\ \ observation of a photon as well as the observation of no photons reveals information 
about the system, thereby resulting in a so-called environment-induced measurement 
PQ. These measurements can assist quantum computational tasks, like the controlled 
generation of entanglement, in many different ways 0, [3]. Measurement-based state 
preparation schemes are in general relatively simple and promise very high fidelities. In 
fact, their performance is widely independent of the concrete size of the experimental 
parameters. In general, it is only limited by the accuracy with which the relevant 
measurement outcomes can be detected. 

For example, no-photon measurements can improve the performance of adiabatic 
passages by stabilising the trajectory of the system from one quantum state into 
another [H [5]. They are also able to restrict the time evolution of a system onto 
a decoherence-free subspace, thereby resulting in very robust and relatively simple 
quantum computing schemes [6j [7J [8j [9j [10]. Another way to take advantage of 
dissipation is to use photon measurements for the build up of highly entangled states 
for distributed quantum computing [TH [121 [131 EH]- The feasibility of this approach 
has already been demonstrated experimentally jT5], [16l [TTJ . Moreover, continuous 
photon measurements can be used to control entanglement via direct quantum feedback 

[HI HSl 1201 EH 122] . 
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Alternatively, Verstraete et al. [23] recently suggested to engineer dissipative 
processes to prepare the ground states of frustration-free Hamiltonians without having to 
register measurement outcomes. These processes can be employed for efficient universal 
quantum computing. At the same time, Kraus et al. |24J showed that dissipation can 
be used to prepare multipartite entangled states with efficient relaxation times in the 
number of qubits by designing the system interactions and environmental couplings 
such that the desired state is the stationary state of the system. In general, this 
requires the design of non-local jump operators. Only for certain target states, jump 
operators are found which need to act only on a few neighboring qubits. The scheme [23] 
moreover requires that the desired stationary state of the system is unique. If a system 
possesses more than one stationary state, small perturbations might result in sudden 
jumps between them. Examples are systems with macroscopic light and dark periods 
[251 26\. These can be used to prepare highly entangled states upon the detection 
of a macroscopic dark period [HI [201 121] • See Ref. [27] by Ticozzi and Viola for a 
general framework for the characterisation of attractive quantum Markovian dynamics 
and Refs. [281 [211 [301 EI] f° r other related experimental proposals. 

In this paper we discuss how to cool qubits into the ground state of an applied 
interaction Hamiltonian Hi Qt . The application of the proposed state preparation scheme 
to concrete physical systems of interacting atomic qubits, like atoms in optical lattices, 
ion traps, or atoms in optical tweezers, is straightforward, since we do not require 
complex system-reservoir interactions which result in non-local jump operators. Instead 
the cooling process is realised via laser driving of auxiliary atomic states and the free- 
space emission of photons. As we shall see below, our scheme is analog to laser sideband 
cooling. All qubit states other than the ground state of H\ nt are resonantly driven by 
laser fields. When combined with dissipation, the result of this driving is the transfer 
of an arbitrary initial state into a highly entangled pure state [32]. Cooling atoms 
into entangled states promises high fidelities as long as the relevant coupling constants 
between qubits are much larger than the effective spontaneous decay rates of excited 
atomic states. 

There are five sections in this paper. In Section [2] we introduce the qubit system 
and the cooling device considered throughout the paper. In Section [HI we describe the 
basic mechanism which transfers the qubits into the ground state of H\ nt . In Section H] 
we calculate the fidelity of the prepared state as a function of the system parameters 
and determine the corresponding cooling rates for concrete examples with one and two 
qubits. Finally, we summarise our results in Section [51 

2. Theoretical model 

In the following we consider a quantum system consisting of N interacting atomic qubits 
and a laser cooling device. We introduce the notation that will be used throughout the 
paper and derive the Hamiltonian and the master equation of such a system. 
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2.1. The qubits and their interaction 

Suppose N qubits are encoded in the long living electronic ground states \go) = |0) and 
\gi) = |1) of N atoms in the presence of an interaction Hamiltonian H\ nt . The total 
Hamiltonian of the system reads in the Schrodinger picture 

H = -f/Free + -f^Int , (1) 

where H-pree is the free Hamiltonian of the system. Here we denote the energy of the 
states \g ) and \gi) by Jtwq and fno\. Hence 

N 1 

Hpiee = ^2^2 f^j \9j)ii(9j\ ■ (2) 
i=l j=0 

Moreover, the interaction Hamiltonian H\ nt can be written as 

2^-1 

-Hint — "^2 ^ n \^n)(^n\ > (3) 
n=0 

in terms of its eigenvectors and eigenvalues |A n ) and A n . 

A system with this Hamiltonian can be realised, for example, by trapping single 
atoms in the individual sites of an optical lattice. In this case, interactions between 
neighboring atoms can be induced by superposed laser fields which result in state- 
dependent tunneling rates and level shifts. Indeed it has been shown that it is possible 
to generate a variety of two and three-body interactions in this way [33]. The aim of 
this paper is to analyse a potential cooling mechanism to transfer the qubits into the 
ground state |Ao) of Hi Qt . In case of a ground state degeneracy, the system is cooled 
into a mixed state. 



2.2. The cooling device 

Our cooling device consists of a set of laser fields which couple the atomic ground states 
|<7o) and \g\) to the auxiliary excited states |eo) and |ei), respectively, as shown in 
Figure [TJ We denote the energies of |e ) and \e\) by ftw and Kb\. Since we added two 
levels to each atom, the Hilbert space of the system is now of dimension 4^. In the 
presence of the laser driving, the Hamiltonian in Eq. (TjQ) hence becomes 

H = H Flee + Hj nt + -^Laser (4) 

with i^Free n o w given by 

N 1 

Hpree = ^2 ^2 ^ l&M&l + ^3 \ e j)n( e j\ ' ( 5 ) 
i=l j=0 

Notice that the interaction Hamiltonian Hi nt can no longer be written as in Eq. ([3]). In 
addition to the interactions between qubit states, couplings between states with at least 
one atom in the excited state have to be taken into account. 

To do so, we introduce additional states |A n ) with n ranging from 2 N to 4 N — 1. 
Suppose, these states form an orthonormal basis together with the 2^ qubit ground 
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Figure 1. The system consists of a collection of interacting atoms with ground states 
|<7o) and \gi) which encode one qubit. Moreover, K laser fields with different frequencies 
are applied which couple these states to the excited states |eo) and |ei) with the 
shown spontaneous decay rates. For simplicity we assume that the states \gj) and \ej), 
respectively, are of the same energy. Here fij 1 '^ and denote the Rabi frequency 
and the detuning of laser k with respect to the \gj) - \ej) transition in atom i. 



states |A n ) in Eq. ([3]). Multiplying Hi nt with the identities ^^"q 1 |A m ) (A m | and 
Yln^ 1 l-^n)(A n | from the left and from the right, respectively, we then find that it can 
be written as 

2 N -1 4 N -1 4 W -1 

Hint = ^A n |A n )(A n | + 2j Xmn [A m )(A n | ■ (6) 

n=0 m=2 N n=2 N 

For reasons which become obvious later we do not assume that the newly introduced 
states |A n ) with n > 2 N are eigenvectors of Hi Qt . The second term in Eq. ([6]) therefore 
contains non-diagonal terms with the x mn = (A m |ifi nt |A n ) being coupling coefficients 
and diagonal terms with the x nn = (A n |ifi nt |A n ) describing energy shifts. 

Suppose the state |^o) and and the states |eo) and |ei), respectively, are of the 
same energy, i.e. 

ujq = uji = uj and ujq = u)\ = Co . (7) 

In this case, the detunings of the applied laser fields depend neither on i nor j. In 
the following, we denote the Rabi frequency of laser k with respect to the \gj) - |ej) 
transition in atom i by Q^' k \ If K different laser fields are applied, then H^ aser equals 
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in the usual rotating wave approximation 

NIK 

#Lasc r = E E E ^ m t k) e< Q -^ Ak > MM + H.c. (8) 

i=l j=0 k=l 

In general it is not possible to find an interaction picture which removes the time 
dependence from this Hamiltonian. 



2. 3. The effect of spontaneous emission 

In the following, we use the master equation to model spontaneous emission from the 
excited atomic states \ef). If 1^- denotes the corresponding decay rate, then the density 
matrix of the atoms evolves according to 

p= - l -[H,p] + C{p) (9) 

where the Lindblad operator C is given by 

N 1 i i 

^ = Y.Y. T ^p^ -\^^Ip-\p< ^ ]■ ( 10 ) 

i=l j,k=0 

The reset operators Bftf., defined as 

42 = 1^)^1, (ii) 

model transitions from |efc)j into \gf)i with Tj k being the respective spontaneous decay 
rate. The overall decay rate of \ek)i equals Tj = Tj Q + Tji. 

Although the reset operators in Eq. ( ITTj) are local, the state of the atoms 
immediately after an emission can be highly entangled |34j . The reason for this is 
that the state of a single system after a photon emission is a function of the direction 
of the emitted photon. In other words, the state of the atoms immediately after an 
emission can be almost any state. In many cases, the reset state has some overlap with 
the state |A ) which we want to prepare. As we see in the following, it is the interplay 
between the interaction Hamiltonian ifj nt and spontaneous emission which is responsible 
for the eventual cooling of the atoms into |Aq). 



3. The cooling mechanism 

The state preparation scheme proposed in this paper is analogous to the well-established 
technique of laser sideband cooling. Both schemes require non-resonant laser driving 
and the spontaneous emission of photons. In both cases, the target state is one which 
is easy to reach but difficult to leave. After a certain transition time, one can therefore 
be sure that the system is with a very high fidelity in this state. Since the idea of state 
preparation via cooling can be understood easily by comparison, we start this section 
with a short overview on laser sideband cooling. We then refer to the one-qubit case 
and consider a generalisation of the proposed state preparation scheme to N qubits. 
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Figure 2. Level scheme for sideband cooling of a single atom which is based on the 
resonant driving of the |l,m) - |2,m — 1) transitions, while all other transitions are 
detuned by the phonon frequency v. The ground state |1, m = 0) is therefore the only 
off-resonant state. Crucial for the cooling to work is moreover spontaneous emission 
of photons from level 2 which preserves the phonon number m. 



3.1. Laser sideband cooling of a single atom 

Sideband cooling is an experimental technique commonly used to transfer single trapped 
atoms and ions close to absolute zero temperature [351 ES] . This is achieved with the 
help of an appropriately detuned laser field which couples the electronic states |1) and 
1 2) with coupling strength rjQ to their quantised motion. For simplicity we consider only 
a one- dimensional trapping potential and denote its phonon states by \m). If the laser is 
red-detuned and its detuning A equals the phonon frequency u, then atomic transitions 
which reduce m by one are in resonance, as shown in Figured Simultaneously, unwanted 
transitions which increase m or keep it constant are out of resonance. Suppose, the atom 
is initially in its ground state and has exactly m phonons. The laser then couples this 
|l,m) state to \2,m — 1) and \2,m) with the transition into \2,m — 1) being the most 
effective. If the atom now emits a photon, its state changes into \l,m — 1) with the 
population in |1, m) being almost negligible. Compared to the initial state, one phonon 
is permanently lost. 

Efficient sideband cooling requires that the laser detuning, i.e. the phonon frequency 
v, is much larger than the laser Rabi frequency Q and the decay rate F of level 2, 

v > Q and T . (12) 

In this case, transitions out of the zero-phonon state |1,0) which is the only non- 



Cooling atoms into entangled states 



7 



resonantly driven state are very unlikely. Under ideal conditions, this is the final state 
prepared in the cooling process. Experiments show that its fidelity can be well above 
99 % [35j EE] . Conditions analog to Eq. ( jl~2l) are in the following imposed on the proposed 
scheme for cooling atoms into an entangled state. 



3.2. 1-qubit case 



In the case of only a single qubit, the preparation of the ground state |Ao) of H\ nt requires 
only a single laser field with detuning A. The relevant four-level scheme is shown in 
Figure [U For simplicity we assume in the following that the (real) laser Rabi frequencies 
Slj and the spontaneous decay rates Tj k are the same for both transitions, i.e. 



e n and T 



Ir. 

2 



(13) 



Taking this into account and transferring the Hamiltonian (jlj) into the interaction picture 
with respect to H = Hp iee — ftA(|eo)(eo| + |ei)(ei|) we obtain the time-independent 
Hamiltonian 



1 i 

Hi= ^2 7j m \ e j)(9j\ + H.c. + hA\ej)(ej \ + H lnt 



(14) 



3=0 



To transfer this Hamiltonian into a more useful form, we introduce the excited atomic 
states | A2) and | A3) as 



I An+2) 



' 1 

E 

.3=0 



I An) 



(15) 



Up to normalisation, these are the states i^Lascr | Ao) and -f^Laser |Ai) which couple to the 
eigenstates |A ) and |Ai) of Hi nt via laser driving. Notice that the states |A 2 ) and |A 3 ) 
are of the same energy as long as H\ nt has no effect on the excited atomic states |eo) 
and |ei). Since also the ground states |go) and \g\) are of the same energy (cf. Eq. ([7])), 
the Hamiltonian (fT4"|) of the system in the interaction picture can be written in terms 
of the |A n ) states as 



Hi= \ m l A n)(A n+2 | + H.c. + h{X n - A) |A n )(A r 



(16) 



ra=0,l 



up to an overall level shift hA with no physical consequences. Moreover, the Lindblad 
operator in Eq. (jTOl) now becomes 

- R) k Rj k p - - 



x ^ x ^ 1 

2 



2. 2.- r 



j=0,l k=2,3 

with the new reset operators 



R jk P Rjk ~ 9 R jk R 3k P ~ TsP R )k R 3k 



(17) 



R 



I A,) (A, 



In the derivation of this equations we took advantage of Eq. ([TBI which assumes equal 
spontaneous decay rates for all transitions. 



Cooling atoms into entangled states 



8 




Figure 3. Level scheme for cooling of a single atom into the eigenstate |Ao) of H\ nt . 
The laser drives only the |Ao) - | A2) transition with detuning A\ and the |Ai) - | A3) 
transitions with zero detuning. The excited states | A2} and | A3) can decay into | A2) 
and I A 3 ). 



Suppose the frequency of the applied laser field equals Co — uo — Ai which implies 
A = Ai. Then the ground state |Ai) is resonantly driven, while |Ao) experiences the 
detuning 

A A = A -Ai, (19) 

as shown in Fig. [3J If A^ is large compared to Q and T, then the comparison with laser 
sideband cooling (cf. Eq. (fl2l) ) suggests that the system reaches |A ) after a certain 
transition time with a very high fidelity. Notice that the detuning A^ of |Ao) which we 
need to prepare the target state comes exactly from the fact that |Ao) is the ground 
state of the system. 

3.3. Generalisation to N qubits 

Let us now have a closer look at the case of N laser-driven atomic qubits. For N atoms, 
the relevant state space is of dimension 4^. Our task now consists of finding laser fields 
which couple the eigenstates |A n ) of Hi Qt with n between 1 and 2^ — 1 resonantly to 
excited atomic states while |Ao) remains off-resonance. Achieving this might require up 
to 2^ — 1 laser fields since there are 2 N atomic ground states. Choosing the right laser 
frequencies requires a detailed knowledge of the structure of Hi nt , since this Hamiltonian 
acts also on states with one atom in |eo) or \ei) and causes level shifts and interactions 
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Figure 4. Relevant level scheme for a the cooling of N qubits into the ground state of 
.Hint. The states |A ra ) with < n < 2 N are eigenstates of H\ nt . Moreover, 2^ — 1 laser 
fields couple |A„) with Rabi frequency 2%//i to |A 2 n + „) with one atom in an excited 
state. The different colors indicate different laser frequencies. These are chosen such 
that I Ao) is the only state without resonant driving. For simplicity, the figure does not 
show the off-resonant driving of the other qubit states. All qubit states are possible 
reset states in case of the spontaneous emission of a photon. 



among them. Here we do not need to consider states with more than one atom in an 
excited state as long as these have sufficiently large spontaneous decay rates. As already 
mentioned above, the emission of a photon in general transfers the atoms into a states 
which has some overlap with |Ao). 

In the iV-qubit case and in the presence of the cooling lasers, the Hamiltonian of 
the system is given by Eq. (TjJ with i^Free as in Eq. ([5]), H\ nt as in Eq. (El), and i^Laser 
as in Eq. ([8]). As in the one-qubit case, it is useful to express this Hamiltonian as a 
function of the |A n )-states introduced in Section [2721 In analogy to Eq. ([6j), H can be 
written as 

4*-l 4^-1 

h= E n iA w )(A n | + x mn iA m )(A n | (20) 

n=0 71=0 m^n 

with 

E n = (Aft I if I Aft) and Xmn = (A m |if |A n ) . (21) 

The i? n -terms in the Hamiltonian (|20l) are effective level shifts, while the Xmn-terms are 
time-dependent since they correspond to laser-driven transitions. In the case of states 
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with at least one atom in an excited state they also include the effect of the interaction 
Hamiltonian H\ nt . 

As mentioned before, the vectors |A n ) with n between to 2 N — 1 are eigenvectors 
of Hi Qt . From this one can easily see that 

E n = hijj + X n and Xmn = for < n, m < 2 N . (22) 

For simplicity we assume again that the laser Rabi frequencies Q^'^ in Eq. ([8]) are all 
the same and given by Q. In analogy to Eq. (fT5j) . we moreover introduce states |A m ) 
with one atom in |eo) or |ei) and m between 2 N and 2 N+1 — 1 as 

N 1 



|Am) = 7n 



|A m _2Jv)- (23) 



i=0 j=0 

The energy E m of these states is given by 

E m = hw + (A m |-ffint|A m ) (24) 

and they are pairwise orthonnormal. Indeed they complete the basis which contains the 
first 2 N states |A n ) further. From Eq. (jSJ) we see that the Xmn 

with 0<n<2 N <m< 
2 N+1 depend neither on n nor m. For them we can hence assume Xmn = X with 

2^-1 

x (t) = -hVNtte 1 ^-^ . (25) 

k=l 

The reason for this simplification is that the applied laser fields drive all transitions with 
the same Rabi frequency. 

The main transitions involved in the time evolution of the system are shown in 
Fig. HI In order to cool into the ground state |A ) it is sufficient to tune for example 
the frequency of laser k in resonance with the |A&) - \Xk+2 N ) transition. This means, the 
frequency uj^ of laser k should be chosen such that 

fi^k = E 2 N +k — E k . (26) 

This implies HA^ = (A 2 iv +fc |ifi nt |A 2 Jv+A,.) — A&, since A& = uj — uj — uj k - Moreover, we 
notice that the detunings of the applied laser fields with respect to |A ) are given by 

E 2 N — Eq — (E 2 N +k — E k ) = (A2iv|-ffl nt |A 2 Jv) — (A 2 iV + fc|-fflnt|A2JV + fc) 

+ A fc -A . (27) 

If k ranges from 1 to 2 N — 1, then the ground state |Ao) is the only state which does not 
experience resonant driving. 

The comparison with laser sideband cooling suggests that population accumulates 
in |Ao) as long as the effective driving experienced by this state is much weaker than the 
driving experienced by the \X k ) states with < k < 2 N . More concretely, in analogy to 
Eq. (11 2p . the system parameters should be chosen such that 

- \E 2 n - E Q - {E 2 N +k - E k )\ > VNQ and T (28) 
n 

for all k. Notice that this condition poses an upper bound on the achievable cooling rate 
for a given interaction Hamiltonian Hi nt . However, it should be possible to speed up the 
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cooling process and to obtain nevertheless relatively high fidelities by slowly decreasing 
the Rabi frequency Q in time. While Q becomes smaller, the fidelity of the final state 
can become arbitrarily close to unity. Once populated, we expect that the atoms remain 
much longer in |Aq) than in any other qubit state. 



4. Concrete examples 

In this section, we explicitly calculate the achievable fidelities and the corresponding 
cooling rates for the one-qubit case with an arbitrary interaction Hamiltonian H\ nt . 
Afterwards, we present numerical results for the two-qubit case for the concrete example 
of a spin-spin Heisenberg interaction and describe the cooling of the atoms into a 
maximally entangled state. 



4.1. 1- qubit case 

In the one-qubit case, the fidelity and the cooling rate of the proposed state preparation 
scheme can be obtained easily by calculating the stationary state of the system 
analytically. Using the master equation ([9]) for the time-independent Hamiltonian Hi 
in Eq. (fT4"l) and setting p = 0, we find that the matrix elements of the stationary state 
with respect to the basis {|Aq), |Ai), | A2) , | A3) } are given by 



AT 2 + 4Aj + n 2 
4(2P + A 2 + Q 2 ) 

at 2 + n 2 

4(2P + A 2 + Q 2 ) 

n 2 

P22 — P33 



Poo 
P11 



4(2r 2 + A 2 + Q 2 ) 

-A x n + irn 



P02 - P 20 - 2(2r2 + ^ + Q2) , 

Pl3 " P31 " 2(2P + A 2 + Q 2 ) (29) 



and 



P01 = P03 = Pl2 = P23 = PlO = P30 = P21 = P32 = (30) 

with A\ defined in Eq. (Tl9l . The fidelity F for the preparation of |Ao) is therefore given 
by Poo; since the state preparation is complete once the system reached its stationary 
state. This means 

4r 2 + 3fi 2 

4(A 2 + 2r 2 + Q?) • (31) 

As shown in Fig. El this fidelity can be arbitrarily close to unity. As suggested by 
Eq. (I28p . high fidelities are obtained when Aa is much larger than Q and T. 

The cooling rate 7 C gives an indication for how long it takes to realise the above 
fidelity. To calculate it for the proposed state preparation scheme, we note the 




V r n/r 



(a) (b) 

Figure 5. (a) Fidelity for the preparation of the ground state | Aq) as a function of 
the effective detuning A\ for different Rabi frequencies fl. (b) Cooling rate 7 C as a 
function of fi/T for large detunings A\ obtained from Eq. ([35]) . 



probability flux conservation 

7hPoo = 7c (1 - Poo) • (32) 

Here 7h is the heating rate, i.e. the rate with which a system prepared in |Ao) leaves the 
target state. To obtain 7h we note that leaving |Ao) is only possible via the accumulation 
of a small amount of population in | A2) due to non-resonant laser driving followed by 
the spontaneous emission of a photon with decay rate T/2 into |Ai). Hence 

7h = ^rp 2 2- (33) 

Since the laser driving of the |A ) - |A 2 ) transition is strongly detuned, this population 
is more or less constant in time, even before the system reaches its steady state. 
Approximating p 2 2 by its steady state value in Eq. (|29|) and using Eq. (|32|) . we finally 
obtain the cooling rate 

= m 2 (4A^ + 4r 2 + Q 2 ) 

7c 8(A 2 +2r 2 + fi 2 )(4r 2 + 3fi 2 ) ' 1 ' 

For large detunings Aa, this equation simplifies to 

7c ~ 2(4r 2 + 3fi 2 ) • (35) 

As illustrated in Figure this rate no longer depends on A\. Maximum cooling rates 
mainly require an as large as possible value for Q without violating condition fl28|) . 
Moreover, T should not be much larger than Q. Otherwise, the interaction with the 
environment results in continuous measurements which freeze the atom as predicted by 
the quantum Zeno effect [37] and make it impossible to reach the target state. 

4-2. 2-qubit case 



In this section we discuss the preparation of two qubits in a maximally entangled 
state. As a concrete example, we consider a particular case of the spin-spin Heisenberg 




20 40 60 80 100 120 140 

time x r 



Figure 6. Fidelity F for the preparation of the maximally entangled state of two 
qubits as a function of time t for J = 5 T and for different fi's. The small figures show 
F as a function of the coupling constant J/T for £1 = T and as a function of the Rabi 
frequency Sl/T for J = 5 T. 



Hamiltonian and assume 

where J is a coupling constant and c?j = (erf, of, of 



(36) 



This means, we treat the two 
ground states of each atom as a pseudo-spin described by the Pauli matrices (a x , a y , a z ). 
In terms of its eigenvectors |A n ) and eigenvalues A n , the Hamiltonian (1361) can be written 

as 

3 



H lnt = -3hJ\\ )(\ \+J2 h J\^n)(\n\ 



(37) 



n=l 



with 



Ao) = 


(101) 


Ai) = 


(101) 


A 2 ) = 


|00), 


A 3 ) = 


|H), 



1 10 
110 



(38) 

and A = —3hJ, while Ai = A 2 = A3 = hJ. In the following, we calculate the stationary 
state of the system in the presence of the cooling lasers in order to determine the fidelity 
of the state preparation. 

As in the one-qubit case in Eq. (|T3|) we assume that the laser Rabi frequencies and 
the decay rates are for all transitions the same. In analogy to Eq. (fl5l) we introduce the 
states 



|A n+ 4> = ^= 



i=l,2 j=0,l 



I An) 



(39) 
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with one atom in |eo) or |ei). For simplicity and since it is anyway small, we neglect 
population in the states with both atoms excited. In this case, the time evolution of the 
system remains restricted onto an eight-dimensional subspace for which the vectors |A n ) 
with n between and 7 form a complete basis. In the interaction picture with respect 
to H = Hp TCC — hJJ2i=i2^2j=oi \ e j)a( e j\i the system Hamiltonian including the laser 
driving can now be written as 

3 1 

Hi = 2 h ^ n l A «X A n+4l + H.c. + h(X n - A) \X n )(X n \ . (40) 

n=0 

When the laser detuning A equals HJ, then the states |Ai), |A 2 ) and |A 3 ) experience a 
resonant driving, while |Ao) is off-resonance. In the concrete example considered here, 
the time evolution of the system is analog to the one-qubit case. Under the condition of 
sufficiently small Rabi frequencies fl and decay rates T, it is possible to achieve fidelities 
well above 90%, as shown in Figure El 

5. Conclusions 

This paper discusses how to cool a system of atomic qubits into the ground state |A ) 
of an applied interaction Hamiltonian Hi Qt . Our cooling device consists of laser fields 
which couple the qubit states \go) = |0) and \gi) = 1 1) to auxiliary atomic states |eo) and 
|ei) with non-zero spontaneous decay rates. Laser frequencies should be chosen such 
that the driving of |Ao) is off-resonant, while all other qubit states experience resonant 
laser driving. Once spontaneous emission results in the population of |A ), the system 
remains there to a very good approximation for a wide range of experimental parameters 
and independent of the initial state of the system. For simplicity, we assume degenerate 
qubit states and degenerate excited atomic states and equal laser Rabi frequencies and 
decay rates for all possible atomic transitions. In this way, it is easy to change from one 
coordinate system into another and complex notation has been avoided. However, the 
application of the proposed cooling mechanism to a large variety of atomic systems is 
straightforward. 

Preparing the ground state of N qubits of an arbitrary interaction Hamiltonian 
Hi nt can require adjusting up to 2 N — 1 different laser frequencies. However, as we 
have seen in Section 14.21 this number is significantly lower in case of degeneracies. In 
concrete situations, it is even possible to realise the state preparation with the help of 
optimised laser pulse sequences [39] which naturally contain a wide range of frequencies. 
As in laser sideband cooling, high fidelities require that the detunings seen by |Ao) 
are much larger than the spontaneous decay rate of |eo) and |ei) and the effective 
Rabi frequencies which couple the qubit states to states with one atom in |e ) or \ei) 
(cf. Eq. (128p ). Efficient cooling therefore requires an interaction Hamiltonian Hi nt for 
which the difference between the energy of the |Ao) - \X 2 n) transition and the energy of 
the \X n ) - |A n+2 iv) transitions with n > is sufficiently large. This means, similar to 
adiabatic quantum computation [38], the proposed cooling scheme requires a relatively 
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well distinguished ground state. 

Our discussion is in good agreement with a detailed analysis for examples with one 
and two qubits. In the one-qubit case, it is possible to calculate the achievable cooling 
rates and fidelities analytically by simply deferring them from the stationary state of 
the system. While the cooling rate increases in general with the laser Rabi frequency fl 
(cf. Fig. 0(b)), high fidelities are easier to obtain for relatively small f2's (cf. Fig. [5](a)). 
One should therefore decrease the laser intensity during the cooling process in order to 
obtain high speed and high precision. In the two-qubit case, we consider a Heisenberg 
interaction Hamiltonian H\ nt with three degenerate eigenstates and a non-degenerate 
ground state | Ao) - 

Like most measurement-based state preparation schemes, the cooling process is ex- 
pected to be very robust against parameter fluctuations. Nevertheless, there is no need 
to register certain measurement outcomes. An interesting question is the efficiency of 
the proposed cooling scheme for a large number of qubits. Since the interaction Hamil- 
tonian ifi nt acts not only on the qubit states but also on states with one atom excited, 
the energy gap on the left hand side of Eq. ff28]) decreases rapidly as 1/N for M-local 
Hamiltonians with M <C N. To increase the efficiency of the proposed cooling mecha- 
nism for large numbers of qubits, one could create entanglement for example by initially 
cooling only separate cells of a finite size. This entanglement can then be distributed 
via controlled interactions between neighboring cells. 
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